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Applying the microcanonical definition of entropy to a weakly interacting and self-gravitating 
neutralino gas, we evaluate the change in the local entropy per particle of this gas between the 
freeze out era and present day virialized halo structures. An "entropy consistency" criterion emerges 
by comparing theoretical and empirical estimates of this entropy. We apply this criterion to the 
cases when neutralinos are mostly B-inos and mostly Higgsinos, in conjunction with the usual 
"abundance" criterion requiring that present neutralino relic density complies with 0.2 < fl^i < 0.4 
for h ~ 0.65. The joint application of both criteria reveals that a much better fitting occurs for the 
B-ino than for the Higgsino channels, so that the former seems to be a favored channel along the 
mass range of 155 GeV < m ; i < 230 GeV. These results are consistent with neutralino annihilation 
patterns that emerge from recent theoretical analysis on cosmic ray positron excess data reported 
by the HEAT collaboration. The suggested methodology can be applied to test other annihilation 
channels of the neutralino, as well as other particle candidates of thermal WIMP gas relics. 

PACS numbers: 98.62.Gq, 14.80.Ly, 95.35.+d 



I. INTRODUCTION 

There are strong theoretical arguments favoring light- 
est supersymmctric particles (LSP) as making up the 
relic gas that forms the halos of actual galactic structures. 
Assuming that R parity is conserved and that the LSP is 
stable, it might be an ideal candidate for cold dark matter 
(CDM), provided it is neutral and has no strong interac- 
tions. The most favored scenario 0, 0, 0, El HI El consid- 
ers the LSP to be the lightest neutralino (x?)> a mixture 
of supersymmetric partners of the photon, Z boson and 
neutral Higgs boson . Since neutralinos must have de- 
coupled once they were non-relativistic, it is reasonable 
to assume that they constituted originally a Maxwell- 
Boltzmann (MB) gas in thermal equilibrium with other 
components of the primordial cosmic plasma. In the 
present cosmic era, such a gas is practically collision-less 
and is either virialized in galactic and galactic cluster 
halos, in the process of virialization or still in the linear 
regime for superclusters and structures near the scale of 
homogeneity 0, H, . 

Besides the constraint due to their present abundance 
as main constituents of cosmic dark matter (ri^i ~ 0.3), 
it is still uncertain which type of annihilation cross sec- 
tion characterizes these neutralinos. In this paper we 
present a method that discriminates between different 
cross sections, based on demanding that (besides yield- 
ing the correct abundance) a theoretically estimated en- 
tropy per particle matches an empiric estimate of the 
same entropy, both constructed for actual galactic dark 
halo structures. The application of this "entropy consis- 
tency" criterion is straightforward because entropy is a 
state variable that can be evaluated at equilibrium states, 



irrespective of how enormously complicated the evolution 
between each such state might have been. In this con- 
text, the two fiducial equilibrium states of the neutralino 
gas are (to a good approximation) the decoupling (or 
"freeze out" ) epoque and their present state as a virial- 
ized relic gas. Considering simplified forms of annihila- 
tion cross sections, the joint application of the abundance 
and entropy-consistency criteria favors the neutralinos as 
mainly "B-inos" over neutralinos as mainly "higgsinos" . 
These results are consistent with the theoretical analysis 
of the HEAT experiment |lfl llll [l2j which aims at re- 
lating the observed positron excess in cosmic rays with 
a possible weak interaction between neutralinos and nu- 
cleons in galactic halos. 



The paper is organized as follows. In section 2 we de- 
scribe the thermodynamics of the neutralino gas as it 
decouples. The microcanonical ensemble entropy is ap- 
plied in section 3 to the post-decoupling neutralino gas 
to estimate the change in entropy between freeze out and 
present day conditions, leading to a suitable theoretical 
estimate of the entropy per particle. In section 4 we ob- 
tain an empirical estimate of this entropy based on actual 
halo variables, while in section 5 we examine the conse- 
quences of demanding that these two entropies coincide. 
We summarize our results in section 6. 
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II. THE NEUTRALINO GAS 

Before the freeze out, the neutralinos satisfy the_ ther- 
mal equation of state of a non-relativistic MB gas 



P = m xl n xl [ 1+ 2x)' 



T ' 



(1) 
(2) 



where m^i and n^i are the neutralino mass and number 

ah Ao 

density. Since we will deal exclusively with the lightest 
neutralino, we will omit henceforth the subscript ^i , un- 
derstanding that all usage of the term "neutralino* 1 ' and 
all symbols of physical and observational variables (i.e. 
flo, to, p, n, etc.) will correspond to this specific particle. 
As long as the neutralino gas is in thermal equilibrium, 
we have 



.9 



TO 

/2tF 



-3/2 



exp (-x) , (3) 



where g = 2 is the degeneracy factor of the neutralino 
species. The number density n satisfies the Boltzmann 
equation 00 



3Hn 



(4) 



where H is the Hubble expansion factor and (cr|v|) is the 
annihilation cross section. Since the neutralino is non- 
relativistic as annihilation reactions "freeze out" and it 
decouples from the radiation dominated cosmic plasma, 
we can assume for H and (crlvl) the following forms 



H = 1.66 gl /2 — , 

TO p 

(<7|v|> = a + 6<v 2 ), 



T 2 



(5) 
(6) 



where m p — 1.22 x 10 9 GeV is Planck's mass, <?* = g*(T) 
is the sum of relativistic degrees of freedom, (v 2 ) is the 
thermal averaging of the center of mass velocity (roughly 
v 2 oc 1/x in non-relativistic conditions) and the constants 
a and b are determined by the parameters characterizing 
specific annihilation processes of the neutralino (s-wave 
or p-wave) |2j . The decoupling of the neutralino gas fol- 
lows from the condition 



T = n(a\v\) = H, 



(7) 



leading to the freeze out temperature Tf. Reasonable 
approximated solutions of (7} follow by solving for Xf 
the implicit relation 2j 



In 



0.0764 Trip c (2 + c ) (a + 6 b/xf) m 



(ff»fZf) 



1/2 



(8) 



where g^ — g*(TA and Co ~ 1/2 yields the best fit to the 
numerical solution of J3J and (JJJ). From the asymptotic 



solution of Q we obtain the present abundance of the 
relic neutralino gas 



nnh 2 = 



So TO 

^crit/^ 2 



2.82 x 10 w y„< 



TO 

GeV 



,(9) 



where 



no 
So 



0.264 g 1 ^ 2 to p to 



|a/a;f + 3(6 ■ 



l/Aa)/xf, 



(10) 



and iSq 



4000 cm 3 is the present radiation en- 



Pcrit 



1.05 



tropy density (CMB plus neutrinos) 
10- 5 GeVcm- 3 . 

Since neutralino masses are expected to be in the 
range of tens to hundreds of GeV's and typically we have 
Xf ~ 20 so that Tc > GeV, we can use g^ ~ 80 in 
equations JSJl - (|10fl . Equation JHJ shows how Xf has a 
logarithmic dependence on to, while theoretical consid- 
erations 0, 0, 0, 0, IE 111 related to the minimal supersy- 
metric extensions of the Standard Model (MSSM) yield 
specific forms for a and b that also depend on to. Insert- 
ing into ©-EIJ) the specific forms of a and b for each 
annihilation channel leads to a specific range of to that 
satisfies the "abundance" criterion based on current ob- 
servational constraints that require 0.2 < f^o < 0.4 and 
h w0.65 0. 

Suitable forms for (a can be obtained for all types 
of annihilation reactions 0. If the neutralino is mainly 
pure B-ino, it will mostly annihilate into lepton pairs 
through t-channel exchange of right-handed sleptons. In 
this case the cross section is p-wave dominated and can 
be approximated by © with 0, Q 



0, 



m 2 [1 + mf /to 2 



(11) 



where to; is the mass of the right-handed slepton (to; ~ 
to p|) and cti = g 2 /47r ~ 0.01 is the fine structure 
coupling constant for the U(1)y gauge interaction. If 
the neutralino is Higgsino-like, annihilating into W-boson 
pairs, then the cross section is s-wave dominated and can 
be approximated by © with 0, Q 



0. 



7ra 2 (1 — to^/to 2 ) 3 / 2 



2 to 2 (2 -to 2 /to 2 ) 2 



(12) 



where m w — 80.44 GeV is the mass of the W-boson and 
a-2 = g\l^ — 0.03 is the fine structure coupling constant 
for the S'?7(2)i gauge interaction. 

In the freeze out era the entropy per particle (in units 
of the Boltzmann constant k B ) for the neutralino gas is 
given by 



P + P 
nT 



Xf, 



(13) 
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where we have assumed that chemical potential is negli- 
gible and have used the equation of state (JIJ. From © 
and l |l 3(1 , it is evident that the dependence of sj on m will 
be determined by the specific details of the annihilation 
processes through the forms of a and b. In particular, we 
will use i[TT|) and |T2|l to compute S£ from [f%|)- fT3*|) . 



III. THE MICROCANONICAL ENTROPY 

After the freeze out era, particle numbers are con- 
served and the neutralinos constitute a weakly interact- 
ing and practically collision-less and self-gravitating gas. 
This gas is only gravitationally coupled to all other com- 
ponents of the cosmic fluid. As it expands, it experi- 
ences free streaming and eventually undergoes gravita- 
tional clusterin g fo rming stable bound virialized struc- 
tures 0, 0, . The evolution between a spectrum 
of density perturbations at the freeze out and the final 
virialized structures is extremely complex, involving a 
variety of dissipative effects characterized by collisional 
and collision-less relaxation processes [Ijj [T^, UjJ • How- 
ever, the freeze out and present day virialized structures 
roughly correspond to "initial" and "final" equilibrium 
states of this gas. Therefore, instead of dealing with 
the enormous complexity of the details of the interme- 
diary processes, we will deal only with quantities defined 
in these states with the help of simplifying but general 
physical assumptions. 

While the assumption of an ideal (thermodynami- 
cal) gas complying with Maxwell-Boltzmann statistics 
is justified before the freeze out (because the collision- 
dominated interactions are short-ranged), this type of 
thermodynamical formalism, associated with extensive 
entropy and energy and with a well defined thermody- 
namical limit, cannot be applied to present day collision - 
less neutralinos subject to a long-range gravitational in- 
teraction [THl Hfij. In these conditions, the appropriate 
approach follows from the microcanonical ensemble in 
the "mean field" approximation that yields an entropy 
definition that is well defined for a self-gravitating gas 
in an intermediate scale, between the short range and 
long range regimes of the gravitational potential. This 
intermediate scale can be associated with a region that 
is "sufficiently large as to contain a large number of parti- 
cles but small enough for the gravitational potential to be 
treated as a constant" 15]- Considering the neutralino 
gas in present day virialized halo structures as a diluted, 
non-relativistic (nearly) ideal gas of weakly interacting 
particles, its microcanonical entropy per particle under 
these conditions can be given in terms of the volume of 
phase space |l6[ 



In 



(2mE) 3 / 2 V 



(14) 



where V and E are local average values of volume and 
energy associated with the intermediate scale. For non- 
relativistic velocities v/c <C 1, we have V cx 1/n cx m/p 



and E cx mv 2 /2 cx rajx. In fact, since the microcanoni- 
cal description is more fundamental, definition ((14(1 . eval- 
uated at the freeze out, is consistent with © and 113(1 . 
and so it is also valid immediately after the freeze out 
era (once particle numbers are conserved). Since 1(14(1 is 
valid at both the initial and final states, respectively cor- 
responding to the freeze out (sj , Xf , nj) and to the values 
associated with a suitable halo structure (s <h) , x (h) , n <h) ), 
the change in entropy per particle that follows from l(14() 
between these two states is given by 



• s f 



= In 



f_ ( x f \ 



3/2' 



(15) 



where ((13() can be used to eliminate Sf in terms of x^. 
Considering present day halo structures as roughly spher- 
ical, inhomogeneous and self-gravitating gaseous sys- 
tems, it is safe to assume that near the symmetry cen- 
ter the spacial gradients of all macroscopic quantities are 
negligible [l8L fl9| . Even if one assumes a "cuspy" den- 
sity profile, as predicted by numerical simulations, the 
gravitational potential of the halo can be considered to 
be practically constant within any typical region of up 
to ~ 10 pc 3 within the halo core (a scale that is smaller 
than the resolution scale of these simulations Since 
such a typical region fits reasonably well the conditions 
associated with the intermediate scale of the microcanon- 
ical description, we will consider current halo macro- 
scopic variables as evaluated near the center of the halo: 

,W T (h) «(>) 

In order to obtain a convenient theoretical estimate 
of s^. h) from 1(151) . we need to relate nj with present day 
cosmological parameters like Qq and h. Bearing in mind 
that density perturbations at the freeze out era were very 
small (Snf/rif < 10~ 4 , 0,1a El), the density is prac- 
tically homogeneous and so we can estimate it from the 
conservation of particle numbers: = rig (1 + z^) 3 , and 
of photon entropy: g^fSt = g*oSo (1 + Zf) 3 , valid from 
the freeze out era to the present for the unperturbed ho- 
mogeneous background. Eliminating (1 + Zf) 3 from these 
conservation laws yields 



9 m { 



9*o L J o 



f 



20.46 n 



(16) 



ticmb 
1 o 



= 4.29 x 10 



12 



GeV 



and g* = £*(T CMB ) ~ 3.91 and T CMB = 2.7 K. Since 
for present day conditions na/n 1 ^" 1 = po/p^ and po = 
p CT [ f f2o h 2 , we collect the results from ((1611 and write 
(ffe|) as 



x { + 92.78 + In 



/my 
VGeV/ 



h ^0 Pcrit 

GeV J (x f x^>) 3 / 2 p™ 



x { + 81.31 + In 



\GcVJ 



h 2 fl GeV/cm 3 



{xcx^fl 2 pi 



(17) 
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Therefore, given m and a specific form of (ojv|) associ- 
ated with a and 6, equation l|17l) provides a theoretical 
estimate of the entropy per particle of the neutralino halo 
gas that depends on the initial state given by in © and 
(|13|l . on observable cosmological parameters flo, h and on 
generic state variables associated to the halo structure. 



IV. THEORETICAL VS. EMPIRICAL 
ENTROPIES 



If the neutralino gas in present halo structures would 
strictly satisfy MB statistics, the entropy per particle, 
s<. h) , in terms of p ( ^ } = mn ( ^ ) and x^ — mc 2 /(k B T c (h) ), 
would follow from the well known Sackur-Tetrode en- 
tropy formula [2l| 



h) | = --\ 

■ I MB 2 

94.42 + In 



In 



m 4 c 3 



ft 3 (27TxH 3/2 /^ h> 



/my 
VGcvJ 



1 



3/2 



GeV/cm 3 

Pc 



(18) 



Such a MB gas in equilibrium is equivalent to an isother- 
mal halo if we identify 22] 



(19) 



where a 2 h) is the velocity dispersion (a constant for 
isothermal halos). However, as mentioned before, the 
assumption of MB statistics, in which entropy and en- 
ergy are extensive quantities, does not apply for self- 
gravitating colissionless gases. Hence, an exactly isother- 
mal halo is not a realistic model, not only because of these 
theoretical arguments, but also because its total mass 
diverges and its distribution function allows for infinite 
particle velocities (theoretically accessible in the velocity 
range of the MB distribution). Therefore, this case will 
not be considered any further. 

More realistic halo models follow from "energy trun- 
cated" (ET) distribution functions [H |H M ^ 
that assume a maximal "cut off" velocity (an escape ve- 
locity). Therefore, we can provide a convenient empirical 
estimate of the halo entropy, s^. h) , from the microcanon- 
ical entropy definition l|14|) in terms of phase space vol- 
ume, but restricting this volume to the actual range of 
velocities (i.e. momenta) accessible to the central parti- 
cles, that is up to a maximal escape velocity v e (0). It is 
reasonable to assume a relation of the form 

vl(0) = 2 | $(0)| ~ aafjO), (20) 

where $(r) is the newtonian gravitational potential. 



With the help of (|20|l we have then 



,(h)| 



In 



.17 + In 



/ m \ 
VGeV/ 



7^ 



3/2 GeV/cm 3 



Pc 



(21) 



where we used x^' = c 2 /tx ( 2 h) (0) as in ljl§|) . As expected, 
the scalings of 12 1H are identical to those of l]18p. 

In order to provide a numerical estimate for a, we re- 
mark that by a suitable choice of a core radius and a to- 
tal radius, i? core and i? tot , any realistic halo model can be 
accommodated to a generic galactic halo density profile 
consisting of a region of constant density < r < i? corc , 
joined continuously to an outer region i? corc < r < Rt at 
in which p decays as r -2 and goes to zero for r > R tot . If 
we further assume a velocity dispersion c 2 h) (0) ~ 2t) 2 ot , 
where v lat is the maximal observed rotation velocity, the 
parameter a defined in l|20|) will be given by: 



2 + 4 In 



(22) 



From this estimation, the minimum value of a occurs 
for the combination of large i? coro and small i? tot , with 
extreme values given (for a large spiral galaxy) by 10 
kpc and 100 kpc, respectively. Conversely, the maximum 
value of a corresponds to small i? coro and large R tot , with 
extreme values safely estimated as 1 kpc and 300 kpc, 
respectively. Considering these extreme values, we obtain 
the uncertainty range 



11.2 < 



< 24.; 



(23) 



Values of a outside these bounds are unrealistic [26, 27, 
l2Sj . an assertion that can be also based on theoretical 
studies of dynamical and thermodynamical stability asso- 
ciated with ET distribution functions [2II2I I25II29I I30I 
I3H l3^ | and on observational data for normal spiral and 
LSB galaxies as well as galaxy clusters fl^l3al34l35ll3q| . 
In any case, it must be noted that a is an integral feature 
of the halo, highly insensitive to the differential details 
of the density profile. 



V. TESTING THE ENTROPY CONSISTENT 
CRITERION 

The expressions for s|. h) obtained from l|17fl and l|21l) are 
not very useful, since they depend on dynamic halo vari- 
ables, such as p<. h) and x 1 -^ , whose values depend on spe- 
cific galactic systems and are subject to possibly large un- 
certainties. However, since both entropy estimates <|17[) 
and (|21|l must correspond to the same quantity, s^. h) , the 
comparison of these independent estimates leads to the 
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following constraint 



x { = 7.85 



In 



(a Xf) 



3/2 



h 2 n GeV 



(24) 



which does not depend on halo variables other than 
a, hence eliminating the effects of potentially large un- 
certainties in the estimation of these quantities. In 
fact, the constraint (|24|) can be interpreted as the con- 
straint on sc = 5/2 + Xf that follows from the condi- 
tion 4 h) L = 4 h) Lm- Since we can use @ and 10 to 
eliminate h 2 ilo, the constraint (|24|l becomes a relation 
involving only x^, m, a, b, a. This constraint is indepen- 
dent of JHJ), which is another (independent) expression 
for Sf = 5/2 + xt, but an expression that follows only 
from the neutralino annihilation processes. Therefore, 
the comparison between Sc h) | th and s£. h) | em , leading to a 
comparison of two independent expressions for Sf, is not 
trivial but leads to an "entropy consistency" criterion 
that can be tested on suitable desired values of m, a, b, a. 
This implies that a given dark matter particle candidate, 
characterized by m and by specific annihilation channels 
given by 2f through JSJ , will pass or fail to pass this con- 
sistency test independently of the details one assumes 
regarding the present day dark halo structure. This is 
so, whether we conduct the consistency test by compar- 
ing © and J23J or {T7| and (|2J . Since the matching of 
either © and J2U or (0 and (|2U shows a weak loga- 
rithmic dependence on m, the fulfillment of the "entropy 
consistency" criterion identifies a specific mass range for 
each dark matter particle. This allows us to discrimi- 
nate, in favor or against, suggested dark matter particle 
candidates and/or annihilation channels by verifying if 
the standard abundance criterion @ is simultaneously 
satisfied for this range of masses. 
Since we can write l|24f) as: 



ln(/i 2 fJ ) = 7.85-x f + ln (ai f ) 



\3/2, 



(25) 



this constraint becomes a new estimate of the cosmolog- 
ical parameters h 2 flo, given in terms of a structural pa- 
rameter of galactic dark matter halos, a, the mass of the 
neutralino, m, and the temperature of the neutralino gas 
at freeze out, x^. This last quantity depends explicitly 
not only on m, but also on its interaction cross section, 
and hence on the details of its phenomenological physics 
through (JSJ). Notice that in reaching we have used 
(|14|l twice, both in estimating the change in entropy in 
going from the freeze out condition to the present virial- 
ized halo, and in obtaining the actual empirical entropy 
in present day halos. However, in these two instances, 
(|14|l was expressed in terms of two different sets of vari- 
ables, involving Qq and h in the former case (through 
the calculation of l|16l) and H17[) ). and containing only 
halo parameters in the latter. In this way, equating the 
two expressions for the same quantity yields a constraint 
between the two sets of variables: equation l|25|l . 



At this point we consider values for the constants a 
and b that define the interaction cross section of the neu- 
tralino, and use (|25|l to plot Qq as a function of m in 
GeV's. Using h — 0.65 and given the uncertainty range 
of a in l|23|l . we will obtain not a curve, but a region 
in the CIq — m plane. Considering first condition i|12[l . 
corresponding to Higgsino-like neutralinos, leads to the 
shaded region in figure la. On this figure we have also 
plotted the relation which the abundance criterion © 
yields on this same plane. It is evident that within the 
observationally determined range of f2o (the horizontal 
dashed lines 0.2-0.4), there is no intersection between 
the shaded region and the abundance criterion curve. 
This implies that both criteria are mutually inconsistent, 
thus the possibility that Higgsino-like neutralinos make 
up both the cosmological dark matter and galactic dark 
matter can be ruled out. 

Repeating the same procedure for mainly B-ino neu- 
tralinos, yields figure lb. In this case, we can see 
that the abundance criterion curve falls well within the 
shaded region defined by the entropy criterion. Although 
we can not improve on the mass estimate provided by the 
abundance criterion alone, the consistency of both crite- 
ria reveals the B-ino neutralino as a viable option for 
both the cosmological and the galactic dark matter. 

As noted above, the results of figures la and lb are to- 
tally insensitive to the values of halo variables, x { ^ and 
p { c ] , used in evaluating (|21|l and l|17f) . Different values of 
these variables (say, for different halo structures) would 
only result in larger or smaller actual values of s*. h) for 
each halo structure (subject to potentially large uncer- 
tainties). Therefore, our results are not only insensitive 
to these uncertainties, but are also scale-invariant. 



VI. CONCLUSIONS 

We have presented a robust consistency criterion that 
can be verified for any annihilation channel of a given 
dark matter candidate proposed as the constituent par- 
ticle of the present galactic dark matter halos. Since 
we require that the empirical estimate Sc h) | om °f present 
dark matter haloes must match the theoretical value 
s^. h) | th , derived from the microcanonical definition and 
from freeze out conditions for the candidate particle, the 
criterion is of a very general applicability, as it is largely 
insensitive to the details of the extremely complicated 
evolution of the neutralino gas from its freeze out era 
(hence, it is also insensitive to the structure formation 
scenario that might be assumed). Further, the details of 
the present day halo structure enter only through an in- 
tegral feature of the dark halos, the central escape veloc- 
ity, thus our results are also insensitive to the fine details 
concerning the central density and the various models 
describing the structure of dark matter halos. A crucial 
feature of this criterion is its direct dependence on the 
physical details (i.e. annihilation channels and mass) of 
any particle candidate. 
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(a) HIGGSINO CHANNEL 
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(b) B-INO CHANNEL 
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Recent theoretical work by E. A. Baltz et al. ^} con- 
firmed that neutralino annihilation in the galactic halo 
can produce enough positrons to make up for the excess 
of cosmic ray positrons experimentally detected by the 
HEAT collaboration [ll], H2 • Baltz et al. concluded that 
for a boost factor B s ~ 30 the neutralinos must be pri- 
marily B-inos with mass around 160 GeV. For a boost 
factor 30 < B s < 100, the gaugino-dominated SUSY 
models complying with all constraints yield neutralino 
masses in the range of 150 GeV < m^i < 400 GeV. On 
the other hand, Higgsino dominated neutralinos are pos- 
sible but only for B s ~ 1000 with masses larger than 2 
TeV. However, our results show this second option to be 
unlikely and are in agreement with the predictions that 
follow from as we obtain roughly the same mass 

range for the B-ino dominated case (see figure lb) and 
the Higgsino channel is shown to be less favored in the 
mass range lower than TeV's. 

We have examined the specific case of the lightest 
neutralino for the mostly B-ino and mostly Higgsino 
channels. The joint application of the "entropy consis- 
tency" and the usual abundance criteria clearly shows 
that the B-ino channel is favored over the Higgsino. 
This result can be helpful in enhancing the study of the 
parameter space of annihilation channels of LSP's in 
MSSM models, as the latter only use equations JSJ and 
(|9"|l- ljl"U)) in order to find out which parameters yield relic 
gas abundances that are compatible with observational 
constraints 0, 0> 0, 0, 0, El However, equations (JHJl and 
© (|10fl by themselves are insufficient to discriminate 
between annihilation channels. A more efficient study 
of the parameter space of MSSM can be achieved by 
the joint usage of the two criteria, for example, by 
considering more general cross section terms (see for 
example Q) than the simplified approximated forms 
(|llfl and (|12l) . This work is currently in progress. 



FIG. 1: Figures (a) and (b) respectively correspond to the 
Higgsino and B-ino channels. The shaded regions display f2o 
vs m from our entropy criterion l)25|l with the solid curve giv- 
ing Slo from the cosmological abundance criterion in all 
cases for h — 0.65. The horizontal dashed lines give current 
estimates of flo = 0.3 ± 0.1. It is evident that only the B-ino 
channels allow for a simultaneous fitting of both the abun- 
dance and the entropy criteria. 
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